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Abstract 

In this paper we study the gravitational collapse in loop quantum gravity. We consider the 
space-time region inside the Schwarzschild black hole event horizon and we divide this region in 
two parts, the first one where the matter (dust matter) is localized and the other (outside) where 
the metric is Kantowski-Sachs type. We calculate the state solving Hamiltonian constraint and 
we obtain a set of three difference equations that give a regular and natural evolution beyond the 
classical singularity point in "r = 0" localized. 

Introduction 

Quantum gravity, the theory that wants reconcile general relativity and quantum mechanics, is one 
of major problem in theoretical physics today. The lesson of general relativity is that also the space- 
time is dynamical, then it is not possible to study the other interaction on a fixed background. The 
background itself is a dynamical field. 

One of more diffuse theory of quantum gravity is the theory called of "loop quantum gravity" 
. This is one of the non perturbative and background independent approach to quantum gravity 
(another non perturbative approach to quantum gravity is called " asymptotic safety quantum gravity" 
0). 

In this paper we will apply ideas suggested by full loop quantum gravity to a minisuperspace 
model where we will impose symmetries on the full metric to obtain a reduced model. It is possible to 
implement the Dirac quantization program following the fundamental ideas of loop quantum gravity. 

Some interesting results in this theory are related to the problem of space-like singularity. In fact 
it was shown in , 0] , and [S] that it is possible to solve the cosmological singularity problem and 
the black hole singularity problem by using the ideas developed in full loop quantum gravity theory. 
In |7] the problem of black hole singularity has been analyzed also under the contest of " asymptotic 
safety quantum gravity" ; in that paper authors showed that non perturbative quantum gravity effects 
give a much less singular Schwarzschild metric. 

In this paper we would like to study the gravitational collapse of a dust sphere inside the event 
horizon 8 . We consider the phase of the gravitational collapse when all (dust) matter has crossed 
the Schwarzschild radius. In this phase the symmetric reduced metric is homogeneous and we can 
solve completely the quantum Einstein equations for the minisuperspace model. In particular we have 
two regions (Region 1 and Region 2), the first one is where the matter is localized (Region 1) with 
a Friedmann- Robertson- Walker type metric and the other one outside the matter (Region 2) where 
the metric is of Kantowski-Sachs type [5] . We summarize the model in ADM variables following the 
paper of Nambu and Sasaki [HJ. After this we will pass to Ashtekar variables. In Ashtekar variables 
we quantize the system following ideas suggested from full loop quantum gravity pQ. The technology 
used in this paper was developed in the preview papers 0, 5H]> [TTj - @ and fH3] . 
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Figure 1: Spatial section for the gravitational collapse inside the event horizon for the particular 
value of the radial coordinate xo = 7r/2. 



The paper is organized as follow. In the first section we briefly recall properties of the gravitational 
collapse in ADM variables inside the event horizon, r < 2MGn, (Hj- As well known, here the temporal 
and spatial (radial) coordinates exchange their role and so it is possible to study an homogeneous 
space-time. In the second section we recall the Ashtekar formulation of general relativity in terms of 
the gauge connection A l a and of the density triad Ef . In particular we recall the form of the Ashtekar 
connection and of the density triad inside the event horizon. For the symmetric reduced connection 
we introduce the holonomy and we define the classical Hamiltonian constraint in terms of holonomies 
and of the volume operator. At the end of this section we define the boundary conditions on the S 2 
separation sphere (between Region 1 and Region 2 of above) in terms of holonomies. In the third 
section we recall the polymer quantization scheme and we introduce the kincmatical Hilbert space for 
the gravitational collapse. In the forth section we study the dynamics and we solve the Hamiltonian 
constraint inside and outside the matter. 



1 Gravitational collapse in ADM variables inside a black hole 
(classical theory) 

In this section we recall the results of [5J in ADM variables. We consider the space-time region inside 
the black hole horizon and we study the collapse of a dust sphere. The metric inside the matter 
(Region 1)(0 x ^ Xo) is described by a closed Friedmann universe homogeneous and isotropic, 

ds\ = -N 2 (t)dt 2 + R 2 (t)(d X 2 + sin 2 X (d0 2 + sin 2 6d(f> 2 )). (1) 

The metric outside the collapsing star (Region 2)( X ^ X < °°) is °f Kantowski-Sachs type 

dsj = -N%{t)dt 2 + a 2 {t)d X 2 + b 2 {t){d0 2 + sin 2 9d(/) 2 ). (2) 

We have assumed the 3-surface is described by the same radial coordinate x a ^ so m t ne exterior of 
the dust matter region. 
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The action in ADM variables using 0} and J2Jl is 
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where Vi = 4ir J Q dxsin V 2 = 47r£ J** 1 d\ = 47r£(xi — Xo) an d K — 8ttGn (the index a, b = 
1, 2, 3). I introduced C for dimensional reason ([£] = L) in Vi, and xi is a cut-off on the space radial 
coordinate. The spatial homogeneity enable us to fix a linear radial cell and restrict all integrations 
to this cell [3]. To simplify notations we restrict the linear radial cell to xi ~ Xo = 1- 
The corresponding Hamiltonian is 
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where JWdust = Vii? 3 pdust is the constant total dust matter and the momentum conjugate to the 
3-metric q a b is 
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To obtain a correlation between Region 1 and Region 2 we have to take into account the momentum 
constraint on the S -sphere junction surface that is in x = Xo localized. When we impose the condition 
n xj u = on the surface in x — Xo we obtain the following unambiguous junction conditions 



p = Pi?, sin Xo _ Pa = 
1 3Vi V 2 



P 2 = R sin xo 



0. 



(7) 



We can express the first of the relations Q in terms of R, R, 6, 6 and we obtain the new constraints 
set 

RR sin 2 xo 
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The gravitational collapse inside the horizon in ADM variables is completely defined by the fol- 
lowing four constraints [Hj 
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Figure 2: Penrose diagram for the gravitational collapse inside the event horizon (Region 1 and region 
2) and outside the event horizon (Region3). 

Solving the constraints equations we obtain the known results for the classical dust matter gravi- 
tational collapse [Hj- 



2 Gravitational collapse in Ashtekar variables 

In this section we study the gravitational collapse in Ashtekar variables |12| . In particular we will 
express the Hamiltonian constraint inside and outside the matter and the constraints P\ and Pi in 
terms of the symmetric reduced Ashtekar connection ^3] > [HI • 



2.1 Ashtekar variables 

In LQG the fundamental variables are the Ashtekar variables: they consist of an SU(2) connection 
A l a and the electric field Ef, where a, b, c, ■ ■ ■ = 1,2,3 are tensorial indices on the spatial section and 
i,j,k, ■ ■ ■ — 1,2,3 are indices in the su(2) algebra. The density weighted triad Ef is related to the 
triad e l a by the relation Ef — ^e ahc eijk e 3 b e k c . The metric is related to the triad by q a b = e a e l^ij- 
Equivalently, 

v/det^) q ab = Ef E) S ij . (10) 
The rest of the relation between the variables (A l a ,Ef) and the ADM variables {q a b,K a b) is given by 

Ai = Ti + 7 K ab E b j ^ (11) 



where 7 is the Immirzi parameter and is the spin connection of the triad, namely the solution of 

a i 1 i -nj k 
-b] + 6 jk 1 [a e b] 



Cartan's equation: <9[ a eli + e* fc T? e?, = 0. 



The action is 



S= — [ dt d 3 x \-2Tr(E a A a ) - NH - N a H a - N^Gil , (12) 
J Js L J 

where N a is the shift vector, N is the lapse function and N z is the Lagrange multiplier for the Gauss 
constraint Gi- We have introduced also the notation Em = E a d a = EfT % d a and Am — A a dx a — 
A l a T l dx a . The functions Ti 1 Ti a and Gi are respectively the Hamiltonian, diffeomorphism and Gauss 
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constraints, given by 

H(E?, 4) = -4 e- 1 Tr (i* 1 ^ E a E^ - 2 e" 1 (1 + 7 2 ) E^E^K^ 
H b (E?,Ai) = E?F> db -(l + 7 2 )K i a G i 

g i (E?,Ai) = d a Et + e t i j AiE%, (13) 

where the curvature field strength is F a b = d a At — dbA a + [A a , Ab] and the determinant is e = det(e^) = 
v /|det(£ [1] )sgn(det(£ [1] ) 0. 

2.2 The Hamiltonian constraint inside the matter 

In this section we recall the Ashtekar variables for an homogeneous and isotropic space-time of topology 
R x S 3 ^31- For this space the Ashtekar 's connection and the densitized triad are 

Ai = cu I a 5 i I , E?=pXf6i, (14) 

where or are the left-invariant one- forms and Xj are the left-invariant vector fields fulfilling oj 1 (Xj) = 
6 j. The Hamiltonian constraint in terms of the variables (c,p) is 

H = H { E HW) +H M = —c(l- c)sgn(p) y/\p\ + H M (15) 

where H)^" W ' is the Euclidean part of the Hamiltonian constraint inside the matter and Hm — M^ust 
is the Hamiltonian constraint for the dust matter. I recall that the metric of the space-time in the 
region where the matter is localized is 

ds 2 = -dt 2 + a 2 (t){d X 2 + sin 2 x(sin 2 6d<j) 2 + d6 2 )). (16) 
The volume operator for the space section inside the matter in Ashtekar variables is 

V (RW > =27r( Xo ™cosxosin X o)bl i =V( X o)H f . (17) 
The classical theory is define by the symplectic structure {c,p} = |13|. 

Holonomies and the Hamiltonian constraint inside the matter. We introduce the holonomies 
for a space-time homogeneous and isotropic in the direction / 

h (RW) = e5oCTl = cos ( C(5o / 2 ) + 2T /S in( C (5o/2), (18) 

and we express the gravitational part of the Hamiltonian constraint in terms of the holonomies 

H m = _ 8 V T,[hf W \f W \f W) ~\f W) '\\f^-\^ 
K 70 V(xo) IJK 

2.3 The Hamiltonian constraint outside the matter 

We recall that outside the matter but inside the horizon we have a Kantowski-Sachs type space- 
time. An homogeneous but anisotropic space-time of spatial section E of topology £ = K x S 2 is 
characterized by an invariant connection 1-form A^ of the form |14| . |15| 

A m = A(t) t 3 dr+{A 1 (t) ri + A 2 (t)r 2 ) dO + (A 1 (t) r 2 - A 2 (t)n ) sin 6 # + r 3 cos 9 dcj>. (20) 
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The Ti are the generators of the SU(2) fundamental representation. They are related to the Pauli <Ji 
matrix by Ti — — \oi- On the other side the dual invariant densitized triad is 



E [x] =E{t)r s sinO^ + {E 1 (t)T 1 + E 2 (t)T2) sinfl — + (E\t) t 2 - E 2 (t)n) — . (21) 



The coordinate r is related to the coordinate x used in the first section by r = IpX- 

We are interested to the Kantowski-Sachs space-time with space section of topology 1 x S 2 ; the 
connection Am is more simple than in H2U|) with A 2 = A\, and in the triad l|21|) we can choose the 
gauge E 2 — E 1 19 . There is a residual gauge freedom on the pair [Ai.E 1 ). This is a discrete 
transformation P : {A\, E 1 ) — > {—A\, —E 1 ); we have to fix this symmetry on the Hilbert space. The 
Gauss constraint is automatically satisfied and the Euclidean part of the Hamiltonian constraint is 



= 8nlpV2sgn{E) r 
1WI 



2AEA 1 E 1 + (2(A 1 ) 2 - l)^ 1 ) 2 . (22) 



We recall that the relation between the metric and density triad formulation and that the volume 
of the space section S are 

q ab = 6dag(^^-,\E\,\E\ sin 2 fl) , V = 4W2Z P /^|£ 1 |. ( 23 ) 

The phase space consists of two canonical pairs A, E and Ai^E 1 and the simplectic structure is 
given by the poisson brackets, {A, E} = j^— and {Ai.E 1 } = j^j— [El The coordinates and the 
momenta have dimensions: [A] — L _1 , [Ai] = L°, [E] = L 2 and [E 1 ] = L. 



Holonomies and Hamiltonian constraint outside the matter. We introduce the fiducial triad 
°ej = diag(l, l,sin _1 9) and co-triad °lo^ = diag(l, l,sin0), and define the holonomy : h — e-> A W = 

e fA [1]a dx« = jAi^ndX = e f A> I °u J I a °e'}u J nd\ = g / A} u 1 T < dA j where u a = dx°_ = (dr^^) and 

u 1 = °u!^ u a . The holonomy along a curve in the direction "7" is given by 

hi = exp J A^Tidx 1 — exp[AfXolp ts], 

h 2 = exp J A^ndx 2 = exp[Ai^ fa +n)], 

h 3 = exp J Aindx 3 = exp[A lf i (r 2 - n)], (24) 

where A\ = (0, 0, A), A\ = (At,At,0) and A\ = (-A 1 ,A 1 , 0). The connection in J23J is integrated in 
the direction "I"; [i Q lp is the length of the curve along the direction r, /x is the length of the curve 
along the directions 6 and <j> @] 1 . The length are defined using the fiducial triad °e a I . Introducing the 
normalized vectors n\ = (0, 0, 1), n 2 — ^(1) 1) 0), n\ = -^.( — 1, 1, 0) we can rewrite the holonomy hj 



1 In the reference |4] the authors obtain (iq ~ 1. 
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in the direction "7" as 

hi = exp (A/ n} n) = cos + 2n) n sin (^y-) , 

hi = cos — - — J + 2r 3 sin — - — J , 

h 2 = cos - J + V2(r 2 + Ti) sm - J , 

ft 3 = cos - J + V2(r 2 - ti) sm - J , (25) 

where Aj =i = Alpfj, and Aj =2 = Ai =3 = Ai^q\/2. 

We now write the hamiltonian constraint ill'L'l) in terms of holonomies 

H E = Yl e " KTr [hihjhfhj 1 ^!^ h^ihx.V}] , (26) 

K7A *0 UK 

where fy/jj = exp(-^ C/j r 3 ) = cos(/i§ Cu/2) -2r 3 sin(^ C/j/2) and C/,/ = 5 2 /5 3 J - ^j^j- 



2.4 Boundary conditions in Ashtekar variables 

In this section we recall relations between the Ashtekar and the ADM variables. In Region 1 where 
the matter is localized the relations are JS| 



(27) 



= 2( 1 " 7 ^)' W=a 



In Region 2 instead the relations can be obtained comparing the Hamiltonian constraint written in 
terms of ADM variables and in terms of Ashtekar variables. From this analysis we obtain 



A ■ a A lb 

A = i — , Ai = —= — 

N 2 V2N 2 



\E\ = b 2 



(E 1 ) 2 



a 2 b 2 



Introducing (jT7|) and in © obtaining 

(2c - 1) sin 2 xo - y f W\A 1 = 0. 



(28) 



(29) 



We can express (|29l) in terms of holonomies using {!%]) and and we obtain 



±- VW\ [4Tr(/i( W) n) + S ] sin 2 ( Xo ) - — VW\^[h 2 (r 2 + n )} + 0( M o) = 0. (30) 
Oo Mo 

Introducing the explicit form of holonomies in (|30|l we obtain the following finally form for the bound- 
ary conditions (JHJ) 



1 



4 sin 



\p\ sin(xo) 



(^°) + So] sin 2 ( Xo ) - ^ 2j2]E\ sin (^|°) + O^o) = 0, 



E = 0. 



(31) 



3 Quantum kinematics 

We want to quantize the collapse inside the horizon using techniques from loop quantum gravity. Now 
we build the kinematical Hilbcrt space TLkin- We define a graph T as a countable number of triple 
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of points (ci,fiEi ) Heh)-> where a, [iEi, Me 1 ; € K< We denote by Cyl r the vector space of functions 
/(c, A, A x ) (/ : R 3 -» C) of the type 

/(c > A > Ai) = 2/« fc e 2 + 2 +^7t^. (32) 

where c, A, Ai € R, Ci,fiEj, ^E 1 k S R , /^fe S C and fc run over a finite number of integers 
(labelling the points of the graph). We call the function /(c, A, Ai) in Cylr cylindrical with respect to 
the graph T. We consider all possible graphs (the points and their number can vary from a graph to 
another) and denote by Cyl the infinite dimensional vector space of functions cylindrical with respect 
to some graph: Cyl = lJ r Cyl r . Thus, any element f(c, A, A\) of Cyl can be expanded as in 

c i[±f£i Ip A _ J 

where the uncountable basis e~ ~ (g> e 2 ® e V2 is now labeled by arbitrary real numbers 
(c,Mb , Mb 1 )- A basis in Cyl is given by Me, Me 1 ) = |m) ® |Me) ® Ime 1 )- Introducing the standard 
bra-ket notation we can define a basis 0] in the Hilbert space via 

(c|m) <g> (A\fj, E ) ® (Ai|^ B i) = e 1 ^ <g> e'"^ P " ® e . (33) 

The basis states l|33|l are normalizable in contrast to the standard quantum mechanical representation 
and they satisfy 

(fj,,fj, E ,IJ, E i\v,v E ,u E i) = <W<W^e <W.«bi' ( 34 ) 

The Hilbert space 7ifei„ is the Cauchy completion of Cyl or more succinctly Tikin = -^(R^o^, d/io), 
where Reo/ir is the Bohr-compactification of R and d/io is the Haar measure on R 3 3o ^ r - 

We can quantize the theory using the standard quantization procedure {,}—> — i [ , ] • We recall 
the fundamental difference from the standard Schrodinger quantization program. In the ordinary 
Schrodinger representation of the Weyl-Heisenberg algebra the classical fields c, A and A\ translate 
in operators. In loop quantum gravity on the contrary the operators c, A and A\ do not exist |17| . 
We can not promote the Poisson brackets to commutators [c,p] = ijlp/3, [A,E] = i-flp/Air and 
[A±,E ] = ijlp/l6ir; rather, the quantum fundamental operators are c, E, E , hi. The momentum 
operators can be represented on the Hilbert space by 

n^-il^— E^-i^— E 1 ^-^— (35) 
P 1 3 dc' 47r dA ' L *167T£L4i- (35) 

It is easy to calculate the spectrum of these two momentum operators on the Hilbert space basis. This 
is given by 

PMiMBiMb 1 ) = 7, Mi ME, Me 1 ), 

6 

c-l \ ^Ellp, , Ai, > HE 1 lip, v , qft N 

E\h,he,Hei) = — 5 IMjMe^e 1 ), E \(M,fj, E ,n E i) = — — 7 =-\fj,,fi E ,ij lE i}. (36) 

o7T 167rv2 

We have to fix the residual gauge freedom on the Hilbert space outside the matter. We consider the 
operator P : /xg, Mb 1 ) ~~ * Im,Me, "Mb 1 ) an d we impose that only the invariant states (under P) are 
in the kinematical Hilbert space. The states in the Hilbert space are : A- [|m, Mb, Me 1 ) + |Mj Mb, — Mb 1 )] 
for fi E i 7^ and the states |m, Mb,0) for = 0. 

The holonomy operators hi in the directions r, 0, <fi of the space section Kx S 2 outside the matter 
are : h!^ E \ %2 E ^ and h^ El \ where /i E lp is the length along the radial direction r and /j El is the 
length along the directions 9 and 4> ( au the length are define using the fiducial triad °ef). Inside the 
matter the holonomy operators are hJj RW ^ where /i is the length along any direction \, 6 or <f> inside 
the matter. The holonomies operators act on the Hilbert space Tikin by multiplication. 
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4 Hamiltonian constraint and quantum dynamics 



In this section we are going to study the dynamics of the model. The general strategy to quantize a 
system with constraints was introduced by Dirac. To implement the Dirac strategy we must impose 
the classical constraints at the quantum level to obtain the physical space. This strategy define the 
quantum Einstein's equations that in general are define by (for 7 = i) 

H\ip) = [-Ae- 1 Tr(F ab E a E b )+H M ] W = 0, 
H b \i/j) = [E]F^ + H Mb \ = 0, 

= [d a E? + e% AiEl + g Mi ] W = 0. (37) 

In the case of our minisuperspace (homogeneous) model the second and the third constraints are 
identically zero then to obtain the physical space we have to solve the first scalar constraint. In 
this paper we want to quantize the model following the full loop quantum gravity ideas. To this 
aim we have expressed the Hamiltonian constraint in terms of homogeneous holonomies and we have 
introduced the polymer representation of the Weyl algebra. At this point we have all the ingredients 
to solve the Hamiltonian constraint inside and outside the region where the matter is localized and to 
impose at quantum level the boundary condition introduced in section T2 .41 

As in non-trivially constrained systems, we expect that the physical states are not normalizable 
in the kinematical Hilbert space. However, as in the full loop quantum gravity theory, we again have 
the triplet 

Cyl c H kin C Cyl* (38) 

of spaces and the physical states will be in Cyl*, which is the algebraic dual of Cyl [Q. A generic 
element of this space is 



The quantum version of the Hamiltonian constraint outside the matter can be obtained promoting 
the classical holonomies to operators and the poisson brackets to the commutators. Using the relations 
in (I25[) we can express the quantum version of the Kantowski-Sachs Hamiltonian constraint (|26[l as 



167ri 



4sin(2a;) sin(2y) ( sm(y)V cos(y) — cos(y)V sm(y) 



+2 (cos(<5)sin 2 (2?/) - sin(<5) (sin 2 (2y) + 2cos(2y))) (sm(x)V cos(x) - cos(x)V sin(z)) , (40) 

where we have introduced the following notations: x — A/iqIp/2, y = ^/2A\iiq/2 and 5 — fJ^/^- 

Using the exponential form for the trigonometric function we can calculate the action of the 
Hamiltonian constraint on the Hilbert space basis (|33J) 6 . At this point we can solve the Hamiltonian 
constraint outside the matter to obtain a first relation for the coefficients ^ #E!(1e1 in the l|39|l . The 

constraint equation He\i/j) = is now interpreted as an equation in the dual space (ip\H E ; from this 
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equation we can derive a relation for the coefficients Vv,me-p b i 



"I"(^A 1 e+2aio,A»e1 — Mo — ^ME+2/io,ME 1— ^Mo) 1 / , A I :A I -B+2/iO:P B i — 2^o 
+ ( 1 ' / AtE-2/iO:P B l +3mo ~~ ^ / Me-2^ ,A 1 e i+Mo) ^A I .A I E-2/i ,A I E l+2Mo 

e+2a'0iA'e 1 + 2 Mo 

sin(//g/2) - cos(^§/2) 



(^ME+MO^'e 1 - 4 M0 _ ^ME— MO.Me 1 ~ 4 Mo) VT*>J*.E)/*®1 ~ 4 M0 
^ V^ME+MOiMe 1 _ ^i"E— MOjMe 1 ) ^f^E>fEl ^ (^ME+MOifE 1 + 4 M0 _ ^ / A'E-A l O,PEl+ 4 A t o) V'p.ME,MEl+ 4 Mf) 

2sin(^i /2) (V^ E + Ml ^ ._2p — V/iE— /io,/i B 1_2 Mo) V'/*,/tE,/iEi _ 2/io 



= 0, 



(41) 



where the volume eigenvalues are defined by 



4V8tt 



Mb I 1Mb 1 



(42) 



The other constraint that we must to impose on the state l|39|l is the Hamiltonian constraint inside 
the matter. If we introduce the holonomies l|18fl in l|19l) we obtain the following trigonometric form of 
the operator 



H 



(RW) 



6(1 — i) ,- <=mq(mq+4) 
e 2 



j + 2(1 + i) e <CW) + 2ie 2iCW5 + 2(1 + i) e 31 ^ -ie 4lcWJ 



V(xo)Jf.7^« 

_|_ e ic Mo _ 2(1 + i) e ic Mo(Mo + l) _ 2 e ic/j (Mo+2) _ + i) e ic MoOo+3) _|_ gic (M0+ 4 ) 

sin (£|^) v(RW) cos (£|o J _ cog ^ V(RW) ^ ^c|o 



(43) 



(the operator c is not defined but the exponential and trigonometric operators are defined on the 
Hilbert space). We have simplified the notation using the identification So = Mo- 

The action of the operator (|4^|l on the basis ll3*3l) is 



MRW) , 



7^p«Mo 

+2 (1 +i) |a* - Mo(mo + 2),fiE,fJ>Ei) + 2i Im ~ Moj Mb, Mb 1 ) + 2 (1 + i) |/x - Mo(mo - 2), Mb 1 ) 
-i |/i - ^ (mo ~ 4),^ B ,/x B i) + |/i + ^ (mo - 4), Mb 1 ) - 2(1 + + Mo(mo - 2),fi E , fj, B i) 
-2 I// + (j%,fi Ei fi E i) - 2 (1 + i) |/i + /x (Mo + 2), Mb 1 ) + |m~ Mo(mo + 4),/x b ,Mb 1 )]- (44) 



We have all the ingredients to calculate the action of the Hamiltonian constraint inside the matter 
on the state From equation (^|(^ W) + H M ) = ((^|Fm = (^|M dust , V/x ^ and (^|-ff M = 
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for m — we obtain another recursive relation for the coefficients ^V^e-Msi 

— * (^M+Mo(Mo+5) — ^+^0(^0+3) 

+ 2 (1 + i) (^+^(^+3) - V^+^o^o-i)) ^+^0(^0+2),^,^! 

+ 2 (1 + i) (V AI+Mo ( Mo _ 2 ) - ^u+Mo(Mo-3)j iPij,+ho(.Ho-2),IJ,e,H e i 
~i (^+^0(^0-3) — ^i+^ (mo-4)) ?/V+/io(mo-4), fj, E ,fj, El 
+ (^-^□(^(,-5) — V^_ Mo ( M[) _3)) '0 (!1 - (!lo ( jUo — 4),^ E ,^ El 

— 2 (1 + i) (V^_^ (^ _3) — V^_^ ( jUo _i)J ^/ , J u- J u (^ -2),/iB,/i B i 

~2 (V M _ Mo ( AIO _i) — ^1-/^0(^0+1)) ^h-hI^he^ei 

-2 (1 + i) (V^-^^o+i) - T/ u _ At0 (^ 0+ 3)) '^ Ai - At0 (^ +2), A iE,/i E i 

)) 

7«|kV(xo)Mo 



+ (^1-/10(^0+3) ^"-/1 (aio+5)J < / ; /i-/io(A'0+4),/lE,M J ; 



0. 



3(2 + 1) 

where we have introduced the volume eigenvalue inside the matter 

V {RW) \fi, Mb 1 ) = V^m^M^Mb 1 ), 
72 „ s 



(45) 



^-V(xo)(^) 5 |Ml 



(46) 



The other constraints to impose are the boundary conditions in (|31|) . The first of the two conditions 
assumes the following operator form 



/'o 



4 sin 



m {—J 



Mo 



sin(xo) - 2 



— sin 1 

Mo 



(^H> 



(47) 



and the action of the operator 147(1 on the state 1)39(1 impose another condition on the physical state 
for the gravitational collapse 



Im - Mo I 



Ho 4 l E 4 l E 1 



M + Mol „, . / IM 



^h+v. ,hb,h*,i -*Mo A/ 1 ^ 1 %,he,h b i )sin 2 (xo) 



MeI 
2vr 



M,Mb,/-'_b1— Mo 



(48) 



The last constraint to impose is the second of equations (|31|l . When we apply this constraint on 
the candidate physical state 1(3911 we obtain 



E [V¥i sin(xo) - \/\E\ 4>h,»b,h b i IMjMe.Me 1 ) 



E 



sin(xo) 



1MB I 
47T 



(49) 



From equation ((49(1 we obtain that the gravitational collapse wave function Vv,/-<e,Me 1 depends only 
on the two parameters m and Me 1 - For sin 2 (%o) — 3/4-7T we obtain 



(50) 



11 



We conclude this section summarizing the quantum dynamical results. In this section we have 
solved the quantum Einstein equations and we have obtained that the gravitational collapse wave 
function must satisfy the difference recursive equations (|41|l . 1)45(1 . I|48|) and then we must impose in 
the result the constraint (|50|l . An important result in our simplified analysis is that all the coefficients 
in the difference equations 1)41(1 . 145(1 and ((48(1 are regular in [i — and in [Ie = where the classical 
Schwarzschild singularity is localized. We can conclude that we have a regular and natural evolution 
beyond the classical singular point. In the next paragraph we will show the consistency of the difference 
equations obtained in this section. 

4.1 Physical states 

In the previous paragraph we have calculated four difference equations (|4*T|l , lj4"5|l , (|4^|l and lf5U|l that 
the wave function for the gravitational collapse must satisfy. Now we analyze this system introducing 
the boundary conditions in the dynamics. We write the physical state in a tensor product form 
for the Region 1 and Region 2 

At this point to simplify the problem we follow the following steps. 

1. The first difference equation is 1(45(1 . this equation gives a relation for the coefficients ip^ intro- 
duced in 1(51(1 : we recall ((45(1 using the simplification [1q — 1 and introducing integer values for 
the independent variable fi, /i = m [1q = m (m G TL 2 ) 

—t (Vm+6 - V m +i) tfm+5 + 2 (1 + i) (V„i+4 - V m ) (fi m +3 + 2 i (V m+ 2 - V m ) f m +l 
+ 2 (1 + i) (Vm-i - Vm-2) ipm-l - i (V m -2 - K71-3) <Pm-3 + (V m +4 - Vm+2) <Pm+3 
-2 (1 + i) (V m+ 2 ~ V m ) (frn+1 - 2 (V m ~ Vm-2) <pm-l - 2 (1 + i) (V m -2 ~ Vm-i) (p m -3 

+ (V m -4 - V m - 6 ) + M dust Vm = 0, (52) 

this difference equation can be solved for ip ma +5 introducing the following initial conditions: 
¥m -5, iPmo-3, frng-i, Vm +i and ip mo+3 (m € 2Z). The component ip mo that is the even 
component in ((52(1 can be calculated using the constraint 1(48(1 with /j,e — as we will stress at 
the point 3. of this paragraph. 

2. The second equation is ((41(1 with = 1, jUgi = 2n and = 21 with n,leZ 

(V2l-2,2n-3 — V2I-2, 2n-l) 4>2l-2,2n-2 + (^2i+2,2n-l — V2l+2,2n-S) 4> 2Z+2,2n — 2 
+ 0^-2, 2n+3 — V2/-2,2n+l) 02/-2,2n+2 + (V2Z+2,2n+l _ V2/+2, 2n+3 J v2/+2,2n+2 

+ (sin(^/2) - COs(/x|j/2)) (V«+l,2n-4 - Va-l,2n-4) </>2i,2n-4 

— 2 (V2;+l,2n — ^2; — l,2n) 02/, 2n + (V2i+l,2«+4 — V2;-l,2«+4j </>2Z,2ri+4 /2 

-2sin(^io/2) (V r A(JS +^ ,^ El -2/io _ V2/_l,2n-2) <f>2l,2n-2 

+ / 2J+l,2n+2 — ^2!-l,2n+2) </>2;,2n+2 = 0. (53) 



2 It is possible repeat the same analysis for the components with m 
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We can solve this difference equation to obtain <f>2i,2n giving the boundary conditions : 4>-2i ,2n, 

4>2l ,2n and 02Z,2n o -4, 4>2l,2n -2, 4>2l,2n , <t>2l,2n +2 VZ. 

The last equation is (|48|l with sin 2 (xo) = 3/4-7T (this is useful to simplify the notation), in this 
equation we introduce \xe — 21 and figi = 2n as in the previous equation 

(3/47r)(V|m - 1| tp m -i - + 1| ip m +i - i x/M ^m) <fc«,2n 

- \/3|I|/7r^ m (^2i,2n-l - 02J,2n+l) = 0. (54) 

(a) For ( = we obtain a relation to express ip m in terms of (p m +i and ip m -i. In particular 
introducing the initial conditions </? mo _i and <y9 mo +i (that are the same initial condition 
introduced in (|52p) we can obtain ip mo from equation 



\/|mo - 1| 9?mo-l - \/|to + 1| ifma + l - i V\ m o\ ¥m = 0- (55) 

At this point we have the even component ip mo useful to solve at the same time equation 
H52J) and equation l|55(l (sec point 1.). 

We observe that if we introduce m = (and obviously I = 0) in equation (|54|l we obtain 
if— i = tpi, cpo decouples from the difference equation; this means that not all the initial 
boundary conditions <p mo -5, ¥>m -3, fmo-U fmo+i, Pm +3 are independent; there is one 
consistency relation between the five initial conditions. 

(b) For I ^ we insert the solution ip m of equation (|52|l in Q54JI and we obtain a second order 
difference equation for 4>2i,2n in the variable n. Using the constraint H50JI for sin 2 (xo) = 3/47T 
we can evaluate equation (|54|l on m = \ie = 21 to obtain 



(3/4-n-) ( i/|m - 1| tpm-i - y/\m + 1| tp m +i - i \/|mjVmJ 



<P2l,2r, 

m=2l 



- y/3\l\/lTip m (<p2l,2n-l - <j>2l,2n+l) = 0. (56) 

m—2l 

Introducing the initial condition (f>2i,2n -i VZ and (f>2l,2n VZ (this second condition can be 
obtained from equation JHSJOj we can calculate 4>2i,2n +i- This observation implies that in 
the variable "n" , from equation (|53|) we can calculate the even components 4>2i,2n and from 
(|56|l we calculate the odd components 4>2i,2n+i of the wave function . 

We conclude the section summarizing the results about physical states. Using the news notations 
we can write the physical states in the following way 

|®£>«,n(2Z,n|, (57) 

where ip m is calculated using the two difference equations (|52(l and l)55|l . and <f>2i,n can be obtained 
from equations l(53")l and 



Conclusions 

In this paper we have studied the gravitational collapse in Ashtekar variables following the paper [H] . 
We have studied the gravitational collapse when all dust matter has crossed the event horizon. In 
this particular region the metric is homogeneous and we have applied the technology developed in 
" loop quantum cosmology" and in " loop quantum black hole" papers 0], ^Sj and [Ej- We have 
divided the space-time region inside the event horizon in two parts, the firs one where the matter is 
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localized and the other one part outside the mater. The space-time metric inside the matter is of 
Friedmann-Robcrtson- Walker type and the metric outside the matter is of Kantowski-Sachs type with 
spatial topology RxS 2 . 

The quantization procedure is induced by full "loop quantum gravity" . We have introduced homo- 
geneous holonomies and we have expressed the Hamiltoninan constraint in terms of such holonomics 
in all the region inside the event horizon. The main result is that the Hamiltonian constraint gives 
a regular difference equation for the coefficients of the physical states which are defined in the dual 
space of the dense subspace of the kinematical Hilbert space. We can summarize this result recalling 
that the quantum Einstein's equations 

(i>\H E (h,V)=0, 

^\{H RW {h RW , V RW ) + H M ) = 0, (58) 

and the boundary conditions are regular in r = where the classical singularity is localized (h, V and 
h RW , V RW are respectively the holonomy and the volume inside and outside the matter). 

An important consequence of the quantization is that, unlike the classical evolution, the quantum 
evolution does not stop at the classical singularity and the "other side" of the singularity corresponds 
with a new domain where the triad reverses its orientation. In this simply model we have solved the 
quantum Einstein equations obtaining recursive equations for the coefficients VV-me,m.ei of the physical 
states 

(■01= VWi^Ei^A^AiiJ 1 !, (59) 

We can interpret Vv,me,m e i as the wave function of the gravitational collapse, after the matter has 
crossed the event horizon. We have also showed the consistency of the four difference equations. In a 
future paper we will study the difference equations (|41|l , l|45f) and l|48|l to obtain a numerical solution 
of those equations for a particular value of the boundary condition on the wave function. 

We hope that the analysis performed here will shed light on the problem of the "information loss" 
in the process of black hole formation and evaporation. See in particular for a possible physical 
interpretation of the black hole information loss problem. 
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